A Precise Determination of Electroweak Parameters in Neutrino-Nucleon Scattering 
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The NuTeV collaboration has extracted the electroweak parameter sin^ dw from the measurement of the 
ratios of neutral current to charged current u and T> cross-sections. Our value, sin^ 0^'°" shoU) _ q 2277 ± 
0.0013(stat) ± 0.0009(syst), is 3 standard deviations above the standard model prediction. We also present a 
model independent analysis of the same data in terms of neutral-current quark couplings. 



Neutrino-nucleon scattering is one of the most precise 
probes of the weak neutral current. The Lagrangian for weak 
neutral current v-q scattering can be written as 
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where deviations from po = 1 describe non-standard sources 
of SU(2) breaking, and ^ are the chiral quark couphngs. 



For the weak charged current, 
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the neutral current e| and e|j each contain an additional term, 
— Qsin^ 9w, where Q is the quark's electric charge in units 
of e. By measuring ratios of the charged and neutral current 
processes on a hadronic target, one can thus extract sin^ 6w 
and Po- 

In the context of the standard model, this measurement of 
sin^ 6w is comparable in precision to direct measurements 
of Mw- Outside of the standard model, neutrino-nucleon 
scattering provides one of the most precise constraints on the 
weak couplings of light quarks, and tests the validity of elec- 
troweak theory in a range of momentum transfer far from 
Mz. This process is also sensitive to non-standard interac- 
tions, including possible contributions from leptoquark and 
Z' exchange [1]. 

The ratio of neutral current to charged current cross- 
sections for either ;y or V scattering from isoscalar targets of u 
and d quarks can be written as [2] 



ct{^v^ n e-(+)x) 



{al + r^-'hll (2) 



where 
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and g\ ji = (e^ jj)^ + (e^ ^)^. Corrections to Equation 2 
result from the presence of heavy quarks in the sea, the pro- 
duction of heavy quarks in the target, higher order terms in 
the cross-section, and any isovector component of the light 
quarks in the target. In particular, in the case where a final- 
state charm quark is produced from a d or s quark in the nu- 
cleon, there are large uncertainties resulting from the mass 
suppression of the charm quark. This uncertainty has lim- 
ited the precision of previous measurements of electroweak 
parameters in neutrino-nucleon scattering [3-5]. 

To reduce the effect of uncertainties resulting from charm 
production, Paschos and Wolfenstein [6] suggested consider- 
ation of the observable: 
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i?^ is more difficult to measure than R'^ , primarily because 
the neutral current scatterings of v and v yield identical ob- 
served final states which can only be distinguished through a 
priori knowledge of the initial state neutrino. 



METHOD 

High-purity v and v beams were provided by the Sign Se- 
lected Quadrupole Train (SSQT) beamline at the Fermilab 
Tevatron during the 1996-1997 fixed target run. Neutrinos 
were produced from the decay of pions and kaons result- 
ing from interactions of 800 GeV protons in a BeO target. 
Dipole magnets immediately downstream of the proton tar- 
get bent pions and kaons of specified charge in the direction 
of the NuTeV detector, while oppositely charged and neutral 
mesons were stopped in beam dumps. The resulting beam 
was almost pure v or V, depending on the charge of the parent 
mesons. Anti-neutrino interactions comprised 0.03% of the 
neutrino beam events, and neutrino interactions 0.4% of the 
anti-neutrino beam events. In addition, the beams of almost 
pure muon neutrinos contained a small component of electron 
neutrinos (mostly from K^^ decays) which created 1.7% of 
the observed interactions in the neutrino beam and 1.6% in 
the anti-neutrino beam. 

Neutrino interactions were observed in the NuTeV de- 
tector [7], located 1450 m downstream of the proton tar- 
get. The detector consisted of an 18 m long, 690 ton steel- 
scintillator target, followed by an iron-toroid spectrometer. 
The target calorimeter was composed of 168 (3m x 3m x 
5. 1cm) steel plates interspersed with Uquid scintillation coun- 
ters (spaced every two plates) and drift chambers (spaced ev- 
ery four plates). The scintillation counters provided trigger- 
ing information as well as a measurement of the longitudi- 
nal interaction vertex, event length, and energy deposition. 
The mean position of hits in the drift chambers established 
the transverse vertex for the event. The toroid spectrometer, 
used to determine muon charge and momentum, also provided 
a measurement of the muon neutrino flux in charged current 
events. In addition, the detector was calibrated continuously 
through exposure to beams of hadrons, electrons, and muons 
over a wide energy range [7]. 

For inclusion in this analysis, events are required to deposit 
at least 20 GeV of visible energy (-Ecai) in the calorimeter, 
which ensures full efficiency of the trigger, allows an accurate 
vertex determination, and reduces cosmic ray background. 
Events with i^cai > 180 GeV are also removed. Fiducial cri- 
teria restrict the location of the neutrino interaction to the cen- 
tral region of the calorimeter. The chosen fiducial volume en- 
hances interactions that are contained in the calorimeter, and 
minimizes the fraction of events from electron neutrinos or 
non-neutrino sources. After all selections, the resulting data 
sample consists of 1.62 x lO*' v and 0.35 x 10® 77 events with 
a mean visible energy (.Ecai) of 64 GeV and 53 GeV, respec- 
tively. 



In order to extract sin^ 9w, the observed neutrino events 
must be separated into charged current (CC) and neutral cur- 
rent (NC) candidates. Both CC and NC neutrino interactions 
initiate a cascade of hadrons in the target that is registered 
in both the scintillation counters and drift chambers. Muon 
neutrino CC events are distinguished by the presence of a fi- 
nal state muon that typically penetrates beyond the hadronic 
shower and deposits energy in a large number of consecutive 
scintillation counters. NC events usually have no final state 
muon and deposit energy over a range of counters typical of a 
hadronic shower. 

These differing event topologies enable the statistical sep- 
aration of CC and NC neutrino interactions based solely on 
event length. For each event, this length is defined by the 
number of scintillation counters between the interaction ver- 
tex and the last counter consistent with at least single muon 
energy deposition. Events with a "long" length are identified 
as CC candidates, while "short" events are most likely NC in- 
duced. The separation between short and long events is made 
at 16 counters (~ 1.7m of steel) for Eca.1 < 60 GeV, at 17 
counters for 60 < -Ecai < 100 GeV, and otherwise at 18 coun- 
ters. The ratios of short to long events measured in the u and 
V beams are: 

i?^^p = 0.3916 ± 0.0007 and R^^^ = 0.4050 ± 0.0016. (5) 

sin^ 9w can be extracted directly from these measured ratios 
by comparison with a detailed Monte Carlo simulation of the 
experiment. The Monte Carlo must include neutrino fluxes, 
the neutrino cross-sections, and a detailed description of the 
detector response. 

A detailed beam simulation is used to predict the ly and V 
fluxes. In particular, a precise determination of the electron 
neutrino contamination in the beam is essential. The ratios 
Rc^p and i?exp increase in the presence of electron neutrinos 
in the data sample because electron neutrino charged current 
interactions are almost always identified as neutral current in- 
teractions. 

The bulk of the observed electron neutrinos, 93% in the v 
beam and 70% in the 17 beam, result from Kf^ decays. The 
beam simulation can be tuned with high accuracy to describe 
t'e and Fe production from charged kaon decay because the 
contribution is constrained by the observed and 1/^ 
fluxes. Because of the precise alignment of the beamline ele- 
ments and the low acceptance for neutral particles, the largest 
uncertainty in the calculated electron neutrino flux is the 1 .4% 
uncertainty in the Kf^ branching ratio [8]. Other sources of 
electron neutrinos include neutral kaons, charmed hadrons, 
and muon decays, all of which have larger fractional uncer- 
tainties (10-20%). Finally, small uncertainties in the calibra- 
tion of the calorimeter and the muon toroid affect the muon 
and electron neutrino flux measurements. Additional con- 
straints from the data, including direct measurements of 
and Vf. charged current events and measurements of i/^ events 
in the 17^ beam (which also result from charm and neutral kaon 
decay) [9] reduce the electron neutrino uncertainties. At the 



highest energies (E^ >350 for i'^ and >180 for ly^), the 
beam Monte Carlo underpredicts the measured flux and is thus 
not used. 

Neutrino-nucleon deep inelastic scattering processes are 
simulated using a leading order (LO) model for the cross- 
section augmented with longitudinal scattering and higher 
twist terms. The cross-section parameterization incorporates 
LO parton distribution functions (PDFs) from charged current 
data measured obtained with the same target and model as 
used in this experiment [10,11]. These PDFs include an ex- 
ternal constraint on a" /a" [11], and make the standard as- 
sumptions that ^■u^p(x) = d nix), d p{x) = ^Un{x)and 
s{x) = s{x). Small modifications adjust the parton densities 
to produce the inherent up-down quark asymmetry consistent 
with muon scattering [12] and Drell-Yan [13] data. A LO 

analysis of '^i'^ N ^+fj,~X events [14] provides the shape 
and magnitude of the strange sea. Mass suppression from 
charged current charm production is modeled using a LO slow 
rescaUng formalism [15] whose parameters and uncertainties 
come from the same high- statistics sample. A model 

for CC production is chosen to match EMC data [ 16] ; it is as- 
signed a 100% uncertainty. A global analysis [17] provides 
a parameterization of the longitudinal structure function, Kl, 
which is allowed to vary within its experimental and theoret- 
ical uncertainties. QED and electroweak radiative corrections 
to the scattering cross-section are appUed using code supplied 
by Bardin [18] and from V6.34 of ZHTTER [19], and un- 
certainties are estimated by varying the parameters in these 
corrections. 

The Monte Carlo must also accurately simulate the re- 
sponse of the detector to the products of neutrino interactions 
in the target. The critical parameters that must be modeled are 
the calorimeter response to muons, the measurement of the 
position of the neutrino interactions, and the range of hadronic 
showers in the calorimeter. Precise determination of these ef- 
fects is made through extensive use of both neutrino and cal- 
ibration beam data. Measured detector parameters are then 
varied within their uncertainties to estimate systematic errors. 

An important test of the simulation is its ability to predict 
the length distribution of events. Figure 1 shows event length 
distributions in the final data sample compared to the Monte 
Carlo prediction for our measured value of sin^ 6w Events 
reaching the toroid, which comprise about 80% of the CC 
sample, have been left out for clarity, but are included in the 
normalization of the data. ExceUent agreement within uncer- 
tainties is observed in the overlap region of long NC and short 
CC events. 

RESULTS 

Having precisely determined i?exp' ^exp' ™d their pre- 
dicted values as a function of electroweak parameters sin^ 0^ 
and pq, we proceed to extract the best values of sin^ and 
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FIG. 1. Comparison of v and v event length distributions in data 
and Monte Carlo (MC). The MC prediction for CC events is shown 
separately. Insets show data/MC ratio comparisons in the region of 
the length cut with bands to indicate the \a systematic uncertainty in 
this ratio. 
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TABLE I. Uncertainties for both the single parameter sin^ Qw fit 
and for the comparison of R" and W with model predictions. 



ments of Mw is 80.45 ± 0.04 GeV [20]. 

For the simultaneous fit to sin^ Qy/ and po. we obtain: 



po- This is done by means of a fit that also includes the slow- 
rescaling mass for charm production (mc) with its a priori 
constraint from ii~ data [14]. i?" is much less sensitive to 
sin^ Bw than W , but both are sensitive to mc and po- 

When fitting with the assumption pQ = \, sin^ 6'vi/ is simul- 
taneously fit with the slow-rescaling parameter mc. Like an 
explicit calculation of this procedure reduces uncertain- 
ties related to sea quark scattering as well as many experimen- 
tal systematics common to both v and v samples. Statistical 
and systematic uncertainties in the sin^ Q-^ fit and in the com- 
parison of W and W with the Monte Carlo prediction are 
shown in Table I. 

The single parameter fit for sin^ dw measures: 



po = 0.99789 ± 0.00405, sin^ 6'vi/ = 0.22647 ± 0.00311, 

(7) 

with a correlation coefficient of 0.850 between the two pa- 
rameters. This suggests one but not both of sin^ (f^ ^^cicW) 
or po may be consistent with expectations. We have also 
performed a two-parameter fit in terms of the isoscalar 
combinations^ of effective^ neutral-current quark couplings 
katn? = + {df^R? at « _20 GeV^, which 

yields: 

{gff = 0.30005 ± 0.00137, {gf f ^ 0.03076 ± 0.00110, 

(8) 
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Leading terms in the one-loop electroweak radiative correc- 
tions [18] produce the small residual dependence of our re- 
sult on Mfop and Mniggs- The prediction from the stan- 
dard model with parameters determined by a fit to other elec- 
troweak measurements is 0.2227 ± 0.0004 [20,21], approx- 
imately 3(7 from our result. In the on-shell scheme, where 
sin^ 9w = 1 — M^/M"^, and where M^y and Mz are 
the physical gauge boson masses, our result implies Mw — 
80.14 ± 0.08 GeV. The world-average of the direct measure- 



with a correlation coefficient of —0.017. The predicted values 
from Standard Model parameters corresponding to the elec- 
troweak fit described eariier [20,21] are (gf)"^ = 0.3042 and 
(5?f )2 = 0.0301. 

In conclusion, NuTeV has made precise determinations of 
the electroweak parameters through separate measurements 



^Due to the asymmetry between the strange and charm seas and to 
the slight excess of neutrons in our target, this result is only sensitive 
to isovector combinations at about 3% of the sensitivity of isoscalar 
couplings. 

^Effective couplings are those which describe observed experimen- 
tal rates when the processes described by Eqn. 1 are calculated with- 
out electroweak radiative corrections. 



of R" and . We find a significant disagreement with 
the standard model expectation for sin^ ^^""^^''"^ In a 
model-independent analysis, this result suggests a smaller 
left-handed neutral current coupling to the light quarks than 
expected. 
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